Nonequilibrium Dephasing in an Electronic Mach-Zehnder Interferometer 
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We study nonequilibrium dephasing in an electronic Mach-Zehnder interferometer. We demon- 
strate that the shot noise at the beam splitter of the interferometer generates an ensemble of nonequi- 
librium electron density configurations and that electron interactions induce configuration-specific 
phase shifts of an interfering electron. The resulting dephasing exhibits two characteristic features, 
a lobe pattern in the visibility and phase jumps of n, in good agreement with experimental data. 
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Introduction. — An electronic analog of optical Mach- 
Zehnder interferometry (E-MZI) has been recently real- 
ized [l|, 0, H, 0| by utilizing edge channels of integer quan- 
tum Hall (IQH) liquids. As it is one of elementary types 
of interferometry, it can serve as an important probe of 
electronic coherence [j| and entanglement [H, 0, Q • 

The E-MZI has a simple setup consisting of two arms 
and two beam splitters. Recent experiments @, H[ on it, 
nevertheless, revealed puzzling behavior that is hard to 
understand within a noninteracting-electron description 
Q; the interference visibility of the differential conduc- 
tance shows bias-dependent lobe patterns, accompanied 
by phase jumps of it at the minima of the lobes. There 
may exist some unnoticed fundamental physics behind it. 

Electron-electron interactions may be important for 
the puzzling nonequilibrium behavior. Interaction effects 
were studied [ljj Ell in the tunneling regime by using 
bosonization methods. In Ref. [Io| . interactions between 
an E-MZI channel and an additional one outside the E- 
MZI were considered and the resulting resonant plasmon 
excitations were proposed as an origin of the puzzle. On 
the other hand, roles of the shot noise of an additional 
detecting channel were addressed [TJ, [l3| to understand 
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similar lobes found in a related experiment 

In this work, we propose an intrinsic mechanism for 
the puzzling behavior, which does not require additional 
channels outside the E-MZI. A key observation is that 
the shot noise at the input beam splitter of the E-MZI 
generates an ensemble of nonequilibrium electron density 
configurations in the two arms. Then the electron in- 
teraction within each arm induces configuration-specific 
phase shifts of an interfering electron, and the ensem- 
ble average of the phase shifts leads to nonequilibrium 
dephasing. The combined effect of the shot noise and 
the interaction results in lobe patterns andphase jumps, 
which agree with experimental data 0, 0|. We use a 
wave-packet picture to describe the nonequilibrium den- 
sity and treat the interaction phcnomenologically at zero 
temperature. The inter-arm interaction is ignored. 

E-MZI setup — The E-MZI consists of two sources 
i = 1, 2, two drains i = 3,4, two beam splitters j = a, b, 
and two arms I — u, d of length Li [Fig. Ufa)]. The 



source i = 1 is biased by |e|V, while i — 2 is unbiased. 
Transmission probabilities at the splitter j are Tj and 
Rj(= 1 — Tj). Each arm consists of a single edge channel 
of the IQH liquid with the linear energy dispersion char- 
acterized by constant group velocity Vp. We first consider 
the simple case with L u = Ld = L, and later discuss the 
case with L u ^ Ld- There are two time scales, electron 
flight time tfl = L/vf through an arm and time sep- 
aration Ty — 2-Kh/(\e\V) between successive injections 
of nonequilibrium electrons (with single-particle energy 
S [0, | e | V] ) from the source 1. Their ratio gives the av- 
erage number of the nonequilibrium electrons in the two 
arms at a given time, 
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Nonequilibrium density configurations. — At zero tem- 
perature, the nonequilibrium state of the E-MZI can have 
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FIG. 1: (a) Schematic E-MZI setup. It has two beam split- 
ters j — a,b with transmission probabilities Tj and Rj , two 
arms I = u,d, two sources i = 1,2, and two drains i = 3,4. 
Only the source 1 is biased, (b-e) Schematic nonequilibrium 
ensemble of electron density configurations in the two arms 
resulting from the shot noise at the splitter a when the two 
arms have two nonequilibrium electrons. The configurations 
can be described by two packets (dashed lines) propagating 
towards the drains as marked by solid arrows. The probability 
P m of each ensemble element m is marked. 
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the two equivalent forms in the noninteracting limit, 
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where U(t) is the time evolution operator, |0) is the 
Fermi sea in equilibrium, and c^(E) creates an electron 
with the single-particle energy E in the scattering wave 
ip(E) incoming from the source 1. dHX = nW) cre- 
ates an electron in the wave packet 1J, Il5 |. <fr(X) oc 



J Q dEexp(— iEX/hvF)ip(E), centered at position x = 
X with packet width W = vpTy. Here, n is an integer 
and we use the convention that both the arms amount to 
the range (-X/2.L/2). Thus for -L/2 < X < L/2, the 
packet center is located in the arm u or d. The equiv- 
alence between the two forms in Eq. ((2|) can be easily 
verified from the identity {d(nW), d'(n'W)} — <5„, n ' and 
the Pauli exclusion. While the former form is more com- 
monly used, e.g., in the linear response regime, the latter 
may provide more insightful understanding of the ensem- 
ble of nonequilibrium electron density configurations and 
the resulting intrinsic dcphasing mechanism, as shown 
below. Later we use a combination of the two forms. 

As time t passes, 4>{X) evolves to <p(X + vpt). When a 
packet arrives at the input splitter a, it propagates into 
the arm u with the probability R a or into the arm d with 
T a . This process describes the shot noise 14, and 
generates a nonequilibrium ensemble of electron density 
configurations in the arms. The ensemble depends on 
N , i.e., on how many packets have significant weight in 
the arms. For example, for N -C 1, only one packet has 
significant weight and the ensemble has two representa- 
tive elements m = 1,2, in each of which the nonequilib- 
rium density appears in the arm u (d) with probability 
P m =i — Ra (Pm=2 — T a ) at any given time. In this 
case, the nonequilibrium state has the form of |4 , ) n int = 



I*, 



where 
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m=l(2)/nint 



contains the packet state with significant weight in the 
arm u (d) and no weight in d (u). Similarly, when two 
packets have considerable weight [Figs. QJb-e)], the en- 
semble has four representative elements m = 1,2,3,4, 
with probability P m — i? 2 , T„, T a R a , and R a T a , respec- 
tively. When the interaction is turned on, each density 
configuration causes different phase shift to an interfering 
electron, resulting in the nonequilibrium dcphasing. 

Lobe patterns in transmission probability. — We focus 
on a weak-interaction regime where the energy relaxation 
rate is sufficiently smaller than t^ 1 [||. In this regime, 
a single-particle energy Eq (S [0, |e|V]) is a well defined 
quantity, and thus one can consider the interference of the 
scattering plane wave \ip(Ea)). Within the arms, \i/j(Eo)) 
shows the superposition, \ip(E )) — r\u(E )) + t\d(E )), 
of the plane waves \1(Eq)) in the arms I = u,d, where 
|?'| 2 = R a and \t\ 2 = T a . The phase accumulation of 
\1(Eq)) is affected by the interaction between \ip(Eo)) and 



the rest of the nonequilibrium electrons. We call the rest 
as the environment of the electron in \^(Eq)). 

To see how the interaction affects the phase accu- 
mulation, we write the total nonequilibrium state as 

|*> = |*)nint = \i>(Eo)) O |$£ )- We wil1 later turn 011 a 

weak interaction between tp(E ) and $_e , and ignore the 
interaction among the environment electrons, as a weak 
interaction may only slightly modify the nonequilibrium 
densities from the noninteracting case. Both tp(Eo) and 
the environment electrons in $£ are injected from the 
source 1 and end in either the drain 3 or 4. As $s 
will be traced out to obtain the transmission probability 
T(Eq) of ip(E ) to a drain, one has freedom of choosing 
either the wave-packet or the plane-wave description for 
$b . We choose the former, since it is more insightful 
and convenient for obtaining the density configurations 
of $>e m the arms, which is essential to describe the 
interaction. In the former, $s is described by the mov- 
ing train of wave packets constructed by single-particle 
states with energy G [0, \ e\V] except Eq; the exclusion of 
E negligibly modifies the packets from 4>(X). As dis- 
cussed before, due to the shot noise at the splitter j = a, 
\$E ) has the superposition, |$ Bo ) = £) m c m \^ Eo ,m), of 
the multiple-packet states, each corresponding to an ele- 
ment m (with P m = | c m | 2 ) of the nonequilibrium density 
ensemble. As the spatial configurations of the moving 
packets are repeated in time with periodicity ry, time 
ensemble average over [0, ry], in addition to the average 
over the nonequilibrium density ensemble, is necessary 
to obtain T(Eq), We remark that the ensemble averages 
cannot be captured by a mean-field approach. 

We turn on a weak interaction between tp(Eo) and <&e - 
It causes the phase shift of ^ described by an operator, 
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Q(Eo)\. Here, 

Uq is the interaction strength (assumed [1 71 ] to be inde- 
pendent of |e| V), to € [0,Ty] is the time-ensemble index, 
and Ni(t) = J x€l n(x,t)dx acts on |$e ) to measure the 
number of the environment electrons in arm I at time 
t. The time ordering operator is dropped in [/ p h as the 
packet dynamics 4>(X + vpt) makes [Ni{t) , Ny (t')\ a o 
number [131 ] . 

The interference signal of the interfering electron with 
energy Eq can be obtained from the reduced density ma- 
trix Tr cnv [L> ph (t )|^)(*|^ h (to)], which is found to be 
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and 
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The trace Tr env over all the orthog- 
onal multiparticle states {l^cnv.i)} of the environment 
electrons is evaluated using the fact that \&E ,m) an d 
\^E ,m'^m) are "severely" orthogonal to each other, i.e., 

($cnv,i| nt=l n{x k )\$E ,m){®E ,m>\ rife' = l H x k')\® e nv,i) 
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is nonzero only for m = m! where fc mx , k' mx are finite 
positive integers; any finite number of local density 
operators can not transform a packet of finite width 
into another orthogonal packet. The off-diagonal part 
of the reduced density matrix in Eq. describes the 
dephasing of \ip(E )} due to the interaction. 



The factor 



can be further evaluated as 
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where AN(t) = N u (t) ~ N d (t). In Eq. (J), we 
ignore the number fluctuation, (AN(t)AN(t')) m — 
(AN(t)) m (AN(t')) m , or (S(t )S(t Q )) m - (5(t )) 2 m , based 
on the observation that it is not a crucial factor for the 
nonequilibrium dephasing [l8| . 

From Eq. ((3]), the transmission probability of the elec- 
tron in ip(Eo) to the drain i — 3 is obtained as 



T = T(Eq) = %+T 1 Dcos($ b +Vd), 
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where T = T a T b + R a R b , T x = 2^T a T b R a R b , $ B 
is the Aharonov-Bohm phase, Dcos($b + rjo) = 
(K[e i$B £ m ^m(e^ (to) )m]) tQ , (• ■ • >to means the ensemble 
average over time to G [0, ry], D is the nonequilibrium 
dephasing factor, and r\r> is the phase shift of T. Notice 
that T is independent of Eq when L u = Ld and that Eq. 
([5]) reproduces the noninteracting result when Uq = 0. 

We first consider the regime of N <C 1 and Upta [% ^> 1. 
Though this parameter regime is unphysical [17| . it is 
nevertheless illustrative since it allows the analytic eval- 
uation of D and rju. In this regime, the packet 4> may 
be approximated as a square packet of extension W with 
constant density of l/W, and the nonequilibrium den- 
sity ensemble has two representative elements, one packet 
partially in the I = u arm with P m =i = Ra or in d with 
P2 = T a , as discussed before. Then, (5(t )) m in Eq. 
has the constant value of (— l) m S, 
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and one finds, using Arg[- 



2nh h ' 
; arctan(Im ■ /Re-), 



D = ^cos 2 5+{R a -T a ) 2 sm 2 S, 
t)d = Arg[cos(5 - i(R a - T a ) sinS]. 



(6) 
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(8) 



Notice that S is proportional to the bias V and that D(8) 
(thus T) shows lobe patterns with periodic minima of 
value \R a - T a \ at 5 = tt/2, 3tt/2, • • • [Fig. IJa)]. The 
minimum value increases from zero as T a deviates from 
0.5. And, the minima are accompanied by phase jumps. 
At T a = 0.5, t\d shows sharp jumps of 7r at the min- 
ima, while it is zero or 7r otherwise. The jumps become 
smeared as T a deviates from 0.5. All these features are 
related to the which-path information [l9l ] of the nonequi- 
librium electrons. As \e\V increases, the phase shift (<5) 




FIG. 2: Nonequilibrium dephasing factor D(8) (solid curves) 
and phase shift t]d(S) of T (dashed). In (a), we plot Eqs. 
(|7I8|) . while (b) and (c) show the numerical results; see text. 
In (b), iV varies from to 2 as S increases up to 1.87T. 



increases (decreases) by 5 with probability T a (R a ) as 
the environment electrons are in arm d (u). At T a = 0.5, 
the two possible phase shifts in the opposite direction are 
balanced in probability, resulting in the sharp jumps in 
i]d at the lobe minima. When T a ^ 0.5, the balance is 
broken, leading to the smearing of the jump. 

Hereafter, we discuss numerical results for the weak 
interaction regime with Uota/h < it 



17(. A quantita- 



tive calculation needs to take account of a larger number 
of nonequilibrium ensemble elements than Fig. [T] might 
imply, since the packet 4> has a rather slowly decaying 
tail and more importantly for given N, about 2N pack- 
ets appear in the arms during tg,. Thus the minimum 
number of ensemble elements for a given N is 2 2N - in 
the calculation, 2 9 elements are considered in the regime 



of N < 3. Such variety (of (e 



;<5(t )\ 



1) modifies the lobes 



[Figs.[2tb,c)]. Interestingly, the lobes in D and the jumps 
in rjo of the N <C 1 case are maintained, though the lobes 
now acquire a decaying envelope and the minimum po- 
sitions of the lobes are shifted. This robustness can be 
understood as follows. Among the 2 2N elements of the 
nonequilibrium density ensemble, those with probability 
P m = R 2 a N or T 2N have (AN) m = N or -N, generating 
the same phase shift of ±<5 at all to's as in the N <C 1 
case, while in the others (AN) m varies with to within a 
range smaller than S. After the ensemble average over 
to, the lobes themselves survive due to the former, while 
the latter has less contribution to T, giving rise to the 
decaying envelope and the shift of the lobe positions. For 
smaller Uotn/h, larger number of packets are involved (at 
a given S), resulting in more rapidly decaying envelope 
and thus smaller number of visible lobes [Fig. [2fc)]. 

Lobe patterns in dl/dV. — From the zero-temperature 
current I = (\e\/h) J |e|V dE T(E ) = (e 2 /h)VT, one 
evaluates dl/dV = (e 2 /h)T [l + F(V) cos($ B + ^(V))]. 
Here, F(V) gives the visibility of dl/dV, {dI/dV\ max - 
dI/dV\ min )/(dI/dV\ max + dI/dV\ min ), and n F is the 
phase shift of dl/dV; remember S cx V and L u — Ld- 

The lobes in T give rise to similar lobes in dl/dV. 
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FIG. 3: Visibility F (solid curves) and phase shift rjp (dashed) 
of dl/dV, as a function of 8. 

At T a = 0.5, the visibility F(V) shows a lobe pattern 
whose minima reach zero and the phase t]f jumps by ir at 
the minima of the lobes while staying constant in other 
regions. The number of visible lobes in F(V) becomes 
smaller for smaller Uota/h, similar to D [Figs. [3£a,b)]. 
The first zero of F(V) appears even for small Uota/h 
where the second lobe of T almost vanishes [Figs. [He) 
and[3^b)]. As T a deviates from 0.5, the minimum values 
of F(V) increase from zero, similarly to those of D, and 
the jumps of rjp become smeared. In this case the jump 
around the first minimum of F(V) is sharper than that 
around the second [Fig. G3jc)]. 

Discussion. — The lobes and phase jumps in Fig. [3] 
qualitatively agree with the experimental data 0, H| at 
T a = 0.5, where N ~ 0.4 - 1 at V = 8 fxV in the E-MZI 
with L ~ 10 /im. and vf — (2 — 5)10 m/s. The first min- 
imum of F(V) occurs around 6 ~ 3w/8 (where N ~ 0.4) 
in Fig. [31 while around V — 8 /xV in the data. From this, 
we estimate Uq ~ 1.5 — 9 fieV for the experimental E- 
MZIs. The dependence of the number of visible lobes in 
F(V) on Uota/h [ItJ indicates that the number depends 
on the magnetic field, disorder, equilibrium electron den- 
sity, etc, as in the data. Note that the minima of F(V) 
look periodic in V in Fig.[3]Ja), although not in general. 

The case with AL = L u — ^ can be under- 
stood from Eqs. (|5I7I8[) . D and T)r> have the same 
forms as in Eqs. (|7l8p except for the replacements, S — > 
aV(L u + L d )/2 and rj D ->• r) D + E AL/hv F - aVAL/2, 
where a = \e\Uotfi/(2Trh 2 VF). Then, the shift of the lobe- 
minimum positions is governed by (L u + Ld)/2, and for 
Uotfi/h < 7r, rjp becomes an undulating function that no 
longer shows sharp jumps. These modifications are how- 
ever negligible for AL <gi W = 2-kUvf / {\e\V), in agree- 
ment with experimental data 0. 

We suggest that the combined experimental analysis 
of both I and dl/dV may be useful since I can provide 
the direct information of the dephasing factor D. And, 
our formalism may be applicable to possible nonequilib- 
rium dephasing in other electronic interferometries 1211 . 
Finally, careful treatment of the number fluctuation 18( , 
ignored in Eq. ([¥]), may give further understanding, as it 



may modify the decaying envelope and the lobe positions. 

In summary, we have shown that the nonequilibrium 
density ensemble, generated by the shot noise at the in- 
put beam splitter, can cause the nonequilibrium dephas- 
ing in the E-MZI. Our result suggests that the experi- 
mental data in Ref. [2j may only be the first of its kind 
with more nonequilibrium quantum effects waiting for 
their discoveries. 
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by KRF (2006-331-C00118, 2005-070-C00055) and by 
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Note added. — During our manuscript preparation, a 
preprint [2(| addressing shot-noise effects in the E-MZI 
was reported. It considers a regime of the interaction 
strength stronger than ours. 
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